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Abstract. This paper presents error analysis of stabilizer free weak Galerkin finite el- 
ement method (SFWG-FEM) for a second order elliptic equation with low regularity 
solutions. The standard error analysis of SFWG-FEM requires additional regularity on 
solutions, such as H?-regularity for the second-order convergence. However, if the solu- 
tions are in H' with 0 « s < 1, numerical experiments show that the SFWG-FEM 
is also effective and stable with the (1 +s)-order convergence rate, so we develop 
a theoretical analysis for it. We introduce a standard H? finite element approxima- 
tion for the elliptic problem, and then we apply the SFWG-FEM to approach this s- 
mooth approximating finite element solution. Finally, we establish the error analysis for 
SFWG-FEM with low regularity in both discrete H!-norm and standard L?-norm. The 
(P,(T), Pj. 4 (6), [P41 (T)]®) elements with dimensions of space d = 2,3 are employed 
and the numerical examples are tested to confirm the theory. 
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1. Introduction 


The weak Galerkin finite element method (WG-FEM) is a useful numerical method for 
solving partial differential equations effectively. The WG-FEM is naturally derived from the 
standard finite element method (FEM) and the most important idea is to use the general- 
ized functions and their weak derivatives which are defined as generalized distributions. 
The WG method is first introduced in [13,14] by Wang and Ye for the second order elliptic 
equations, and a stabilizer term is added to WG-FEM in order to enforce the connection 
of discontinuous functions across element boundaries [7,8]. Then the WG method finds 
applications in diverse areas including parabolic equations [24, 25, 27, 28], Stokes equa- 
tions [10, 12], Maxwell equations [9, 11], biharmonic equation [26], Cahn-Hilliard-Cook 
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equation [5], stochastic parabolic equations [29,30], eigenvalue problems [22,23], and so 
on. Nevertheless, the stabilizer makes the finite element formulations and programming 
complex. To remove the stabilizer term, a stabilizer free weak Galerkin finite element 
method (SFWG-FEM) is introduced in [16] by Ye and Zhang for the second order ellip- 
tic equations on polytopal meshes. The main idea of SFWG-FEM is raising the degree of 
polynomials for weak gradient computation to increase the connectivity of weak function- 
s. The SFWG-FEM is firstly based on a weak gradient definition similar to the standard 
WG method, and for the element (P,(T), P; (e), [P;(T)]^), it has been demonstrated to be 
reliable if j > k--n— 1, where n represents the number of element's edges/faces [17]. Sub- 
sequently, the requirement j > k-- n — 1 has been loosened in [3]. Various configuration of 
(PCT), P(e), [P;(T)]9) with l, j > k leads to different schemes, cf. [1, 20,21]. A novel def- 
inition of weak gradient is presented in [18] to the (P,(T), Pi. 4(e), [Pia (T) element 
that employs fewer unknowns than in [17] while maintaining the same optimal order of 
convergence. Thereafter, SFWG-FEM is also successfully applied to several equations and 
reduces programming complexity [1—-3, 17-21]. 

In order to reach the optimal convergence order for approximating the second order 
elliptic equations, in many published literature on WG-FEM and SFWG-FEM, the solution 
is usually assumed to have at least H?-smoothness [1, 18]. As a result, it is demonstrat- 
ed that the convergence rate are at least O(h) in H!-norm and @(h”) in L?-norm with 
step-size h, respectively. However, numerical experiments show that we have the O(I) 
convergence in ||-||-norm and @(h!*’) convergence in L?-norm when the exact solution 
has only H!^?-regularity (0 < s < 1). The (1 + s)-order L? convergence analysis on WG- 
FEM is accomplished in [15], and the core concern in [15] is only on analysis but without 
numerical experiments. There is no such theoretical analysis for SFWG-FEM, so in this 
paper, we are devoted to studying the convergence analysis on SFWG-FEM with low reg- 
ularity, and to proving the discrete H!-norm and L?-norm convergence rate to be (I) 
and O(h!^*), respectively. The numerical examples are also tested to confirm the theory. 
Our strategy is divided into two steps. Firstly, we use H?-regular Argyris elements [4] as 
the bases of FEM to approximate the second order elliptic equation whose solution has 
H'*-regularity, and then we get the optimal discrete H! and L? convergence order to be 
of G(I?) and O(h'**), respectively. Secondly, we utilize the SFWG-FEM to approximate the 
H?-regular finite element solution, and then we have the @(h) convergence in ||-||-norm 
and @(h*) convergence in L?-norm, respectively. As a consequence, we arrive at the op- 
timal convergence rates to be (I) in |||-||-norm and @(h'**) in L?-norm, respectively, if 
the exact solution of the second order elliptic equation is only H!*-regular, which is an 
important supplementary for the new stabilizer free weak Galerkin finite element method 
theory. 

The paper is organized as follows. In section 2, we apply the Argyris finite element 
approximation to a second order elliptic equation with H!** regularity, and provide an 
optimal order error estimate. In section 3, we introduce an SFWG-FEM and construc- 
t the corresponding scheme. In section 4, we perform an error analysis and prove the 
optimal convergence order for the SFWG-FEM. In section 5, we present some numerical 
experiments to validate the theoretical analysis. In section 6, we summarize the main 
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conclusions. 


2. The variational form and H? finite element approximation 


In this section, we introduce a second order elliptic equation with lower regularity, and 
then study its variational equation and an H? finite element approximation. 

Let Q C Rİ (d = 2,3) be a convex bounded polygonal domain with Lipschitz contin- 
uous boundary 0. Denote by L?(Q) and L?(@Q) the square integrable function spaces 
with inner products (-,-) and (-,-)a9, respectively. For m > 1, denote by H™ = H™(Q) the 
standard Sobolev space with norm ||- ||,, and by ||- || the L?-norm if no confusion occurs. 
Define by H' a subspace of H™ with vanishing boundary values on 02. 

Consider a second order elliptic equation 


—V-(A(x)Vu(x)) = f(x), ind, 


2.1 
u=0, on OQ. On 


We assume 
(H1) f e LQ). 


(H2) A— A(x) = (aij())asa € [W 99 (Q)]**4 is a symmetric matrix-valued function in Q 
and satisfies 
allğll? < ETAGE < ool", VEER, xen, (2.2) 


where £T is the transpose of a column vector £ and a, > a, > 0 are constants. 


The variational formulation for Eq.(2.1) seeks u € H J such that 
(AVu, Vv)=(f,v), Wve Hx (2.3) 


Under assumptions (H1) and (H2), the variational equation (2.3) possesses an unique 
solution u € Hl. which is called a weak solution of (2.1) [4]. Throughout this paper, we 
will always use u to denote the weak solution of (2.1) when no confusion occurs. 

For the second order elliptic equations, the solutions may not lie in H? if the boundary 
ƏN is piecewise smooth, or the function f and A have lower regularity. It is necessary 
to establish the error estimate for the problems with low regularity solutions, so we also 
assume 


(H3) The weak solution u of (2.1) has H!^* (0 <s <1) regularity with respect to x € Q. 


We denote by J, a shape-regular conforming triangulation of the domain Q. (triangles 
if d = 2, tetrahedra if d = 3) satisfying the regularity requirements stated in [14]. Let &, 
be the set of all edges or faces in Z,, and E = &,\0Q be the set of all interior edges or 
faces. Let hy be the diameter of each element T € %,, and let h = mox hr. We take d = 2 


as an example. 
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Let s? cH 5 N H? be a finite element space which is spanned by the Argyris bases [4], 
thus, all elements in S? are piecewise defined polynomials of degree no more than 5. The 
finite element approximation to (2.1) is to seek ti, € 5? such that 


(AV, Vin) = (f, Vn), VV, € SP, (2.4) 
which is equivalent to seeking iij € s? such that 
=P, (V - (AVii;,)) = Pf, (2.5) 


where P, is the orthogonal projection from L?(Q) to S? . The classical finite element analy- 
sis for the second order elliptic equation gives that [6] 


Lemma 2.1. Assume (H1)-(H3) hold. Let Ŭŭ, € s? and u € H!™*™ be the solution of (2.1) 
and (2.4), respectively. Then the elliptical regularity inequality ||ti,||. < C||f || holds, and for 
0 « s < 1, there exists a constant C > 0 independent of h such that 


l| Vu — Và;l| € CAPS |l; 
llu — üjl| € Ch’ [Ff |l. 


3. Stabilizer free weak Galerkin finite element approximation 


The primary idea behind WG-FEM is to approximate functions and their corresponding 
derivatives with the generalizations of functions and associated weak derivatives. For the 
WG methods with stabilizer term, the weak derivatives are defined in the sense of gen- 
eralized distributions [13, 14]. To remove the stabilizer term, it is essential to increase 
the connectivity of a weak function across element boundary, thus, weak derivatives are 
approximated by polynomials of higher order [16, 17, 20]. 

In this section, we first introduce weak gradient operators and then perform stabilizer 
free weak Galerkin finite element approximation to (2.1). The corresponding L? projec- 
tions which are essential to later analysis are also discussed. 


3.1. Discrete weak gradient operator 


A weak function on T € J, is denoted by a pair v = (vg, vp} with v € L*(T) and v, € 
L?(8 T). For any given integer k > 1, denote by P,(T) and P,(e) the sets of polynomials on 
T € J, and e € & NOT with degree no more than k, respectively. Various configurations of 
(P,CT), P,. (e), [P,(T)]9) lead to different weak Galerkin finite element methods, and the 
(P, (T), P... (e), [B CT lige: elements are introduced in our stabilizer free weak Galerkin 
finite element method. 

For the integer r > 1, define a discrete weak function space W,.(T) by 


W,(T):— {w = {Wo, Wp}: Wo E P,(T), w, E€ P,_,(0T)}, 
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where P,. ,(O T) is a piecewise polynomial set on OT that P,_;(0T) := (w:w|, € P, 1(e), Ve € 
ƏT Nép}. By patching over each element T € %,, we construct the WG finite element space 
V, as follows 

Va := {v = {v0. Vp} : vlr E W,(T), YT € Hh, 


and let Ve refer to the subspace of V, that 
V0 := {v = {vos v] : v € Vis Vala = 0]. 
Define a space of vector-valued polynomials with degree no more than r + 1 by 
My, = {q : qlr € [PCS VT € Hh. 
Next we introduce the definitions of new weak gradient operators. 


Definition 3.1 (cf. Ye. X & Zhang S. [18]). For v eW,(T) on T € J, a local discrete weak 
gradient Vy -+1,r is defined as a linear operator from W,(T ) to [P.41(T)]¢ satisfying 


(Vw,r41,7: q); = (Vvo, q)r + (Qu Cv; qu Yo), q : nN) or > Vq € [P5 CP3]8; 


where n is a unit outward normal vector at OT and Q, is the orthogonal projection from 
L?(8T) to P,_,(@T). 


Definition 3.2 (cf. Ye. X & Zhang S. [18]). A global discrete weak gradient V,, „+1 is defined 
as a linear operator from V, UH! to Ma satisfying 


(Var lr = Vw rt17r lr), We Vn U H!. 


To simplify notations, we will use V,, to denote both weak gradient operators defined 
above when no confusion occurs. For any function v € H!(Q), we note that it could take 
Vo = v and v, = v in the preceding formulation. Then, for smooth functions, we notice 
the weak gradient becomes directly consistent with the strong gradient, that is, for any 
function v € H!(Q), 


(Vw; q) = (Vv, q), Vq € [P4 ( T)]4. 


3.2. SFWG approximation and corresponding L? projections 


We define a bilinear operator a(-,-) and the corresponding H-norm ||-||| for the SFWG 
algorithm of the second order elliptic equation (2.1) with the assumptions (H1)-(H2). 


Definition 3.3 (cf. Ye. X & Zhang S. [18]). For any w,v € V$, the bilinear operator a(-,-) 
is defined as 


a(w,v):— M (AV, w, VY) 75 
TEX, 


and equipped with the norm |||-|| in v? which is defined by |||v||? := a(v, v). 
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Thus, we obtain that the SFWG finite element approximation to (2.1) is seeking uj = 
{Up, Up} € v? such that 


a(up, Vp) > (f, vo), VVh = {vo, vp} € vp. (3.1) 


It follows from (2.2) that the bilinear operator a(-,-) is bounded, and we have the following 
lemma. 


Lemma 3.1. Assume (H2), for any w,v € V,°, there exists a constant p > 0 such that 


la(w, v)| < Bull v l. 


Next we introduce some useful L? projections and lemmas for later analysis. 

Let Q, be the orthogonal projection from L?(8 T) to P, ,(8 T) as we defined before, 
and Qs be the orthogonal projection from L?(T) to P,(T). Then, define the orthogonal 
projection Q, : H! — V, such that 


(Q.$)lz := {Qo(¢lr),Qa(blar)}, Ve eH!. 


In our case, denote by Q, the L?- orthogonal projection from [L?(T)]? to [P,., ,(T)]4, and 
denote by Q, the orthogonal projection from [L?(T)]* to [P,(T)]4 for integer k > 0. For 
these orthogonal projections, there hold the following properties. 


Lemma 3.2 (cf. Wang J. & Ye. X Zhang S. [14]). Let J, be a finite element partition of ©. 
satisfying the shape regularity assumptions. Then, for any ¢ € H**1, we have 


Y le -Qobl < pel. 0<k<r, 
TEJ, 


Y (IVS - Qvo?) s co? *llelll,,, 0<k<1, 


TEX, 
where C1, C; > 0 are constants independent of h. 


Lemma 3.3 (cf. Ye. X & Zhang S. [18]). For any $ € H'(T) on the element T € J, there 
holds the following commutative property 


Vo z Q(V¢). 


Remark 3.1. We would like to emphasize that the operator Q,(V¢) is not identical with 
V,,Q4, and we can find that for any $ € H'*! the difference is 


(VwQrd. Dr = (V(Qo$), 4) 7 + (QQ — Qoh),4°N) ar 
= —(Qo$. V -q)r + (Qo$. q: n)oz + (Qo — Qv(QoG),4°2) ar 
—(V$,q)r — (p.q ‘N)ar + (Qoo — Qui(QoQ),q- N) ap + (Qu4.q j N) ar 
=(Q,(V¢),q)— (6 -Q,0,q:n)s, 
+ (Qo - QY(Qo$),q-n)oz, Va € [P (T)]*. 
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For later analysis, we introduce a useful lemma for the projection Q}. 


Lemma 3.4 (cf. Ye. X & Zhang S. [18]). For any w € H"*1, there holds 
llQuw — w|| s Ch’ lwll, 


where C > 0 is a constant independent of h. 


4. Error analysis 


In this section, we will investigate the optimal error estimate for the SFWG approx- 
imation (3.1) under the weak regularity assumption (H1)-(H3). For simplicity we will 
focus on the case where the symmetric positive definite coefficient tensor A in (2.1) is a 
piecewise constant matrix with respect to the finite element partition Z,. The result can be 
extended to variable tensors without any difficulty, provided that the tensor A is piecewise 
sufficiently smooth. 

Throughout this paper, we always use ii, and u, to denote the solutions of (2.4) and 
(3.1), respectively. Then we have e, := u — up € v? that 


En = (u — Ün) + (üy — Qnty) + (Qpitn — uj). (4.1) 


The discrete H! estimates of the first and second terms in (4.1) follow from Lemma 2.1 
and 3.4, respectively. The L? estimates of the first and second terms in (4.1) follow from 
Lemma 2.1 and 3.2, respectively. We will focus on the error analysis of the last terms in 
(4.1). Define an error that 


Cn = 169, €p} :— Qi — Up; (4.2) 
where eg = QoŬp — ug and ey = Qpil, — Up. 
Lemma 4.1. Assume (H1) and (H2) hold. Let ti, € S? c H?N Ht be the standard finite 


element solution of (2.4) and up € V, be the solution of SFWG-FEM in (3.1). Then, the error 
en = Qhŭp — uj in (4.2) satisfies an error equation 


a(n, Vh) = A(Qhity — tin, vy) + la (ins Vn) + lo (lg, Vh), V Vn = (vo, Vp} € VP, (4.3) 
where l,(ü,, vp) and 15(iip, vp) are defined as follows: for any given w € H'*1, define two 


; 0 
functionals on V, 


l(w,v):— ys (A(Vw — Q,(Vw)) -n,QbYo — v») 


Te€4Z, 


Io(w, v) := y» (AVw : n, vo — Qsvo)or; 
TEX, 


aT’ 
(4.4) 


for any v = iv, Vp} € ve. 
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Proof. To derive an error equation, we decompose ep into two parts since the operator 
V ,Qj is not identical with Q,(V $). Apparently, we have a decomposition 


en = (Qyüg — Un) + Up — ug. (4.5) 


Then we focus on the term a(ti;, — uy, vj). 
For any vp = {Vo, Vp} € Vp, from Lemma 3.3 and Definition 3.2, it follows that 


(tin, vy) = >| (AV, Vw vy)r = >, (AQ), Yuva) r 


TEH, TEH, 

= P2 (AQ, (Vin), Vvo) r + z (AQ, (VÜ) : n, Qu(vy — vo))oq (4.6) 
TEH, TEH, 

= x (AV üj, VYo)r — Xi (AQ; (Vü): n, Qvo — Vb) or- 
TEH, TECH, 


For each v, = {vo, vy} € V a(up, vy) = (f, Vo) holds. Notice that we have vo|r = P,(T) C 
Slr for each element T € J, and 1 <r € 5. Testing Eq.(2.5) with vg and using Green's 
formula, we have 


(fv) = >) (Pafs v0) r 


TE, 
=>) (-A(- (vàp).vo),, 
TEX, 
= > (AVüs, VYo)r — 2 (AVü, n, vo — Vb) ar (4.7) 
TEH, TEH, 
= 2. (AVitn, Vvo)z — 2 (AVūn: n,QbYo — Vb) ar 
TEH, TEH, 
= >. (AVitn 1, Vo — QpY0) ar» 
TES, 


where the second last equality follows from the fact that 


Y (Avü, n vp) op = 0. 
TEX, 


Combining (3.1) and (4.4)-(4.7), we obtain the error equation (4.3). 


Lemma 4.2. (Trace Inequality) Let T be an element with its edge e € &. For any function 
V € H'(T), then we have 


Illi? < CF holz + hr ll VIT). (4.8) 
Lemma 4.3. (Inverse Inequality) There exists a constant C = C(r) such that 
IVpllr < CAF IMs, YT € HR, (4.9) 


for any polynomial ^j of degree r on T. 
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Lemma 4.4 (cf. Ye. X & Zhang S. [18]). For any w € H'*! and v = (vg,v,] € V, there 
hold 


Iw, v)| S Ch" |lw]l-+ llv, 
Il5Qw, v)| S Ch Iwll- Mv, 


where C > 0 is a constant independent of h. 


Applying the error equation (4.3), we derive the optimal order discrete H!-norm error 
estimate of e, in (4.2). 


Theorem 4.1. Assume that (H1) and (H2) hold. Let it, € S? c H?N H? be the standard 
finite element solution of (2.4) and u, € V, be the stabilizer free weak Galerkin finite element 
solution of (3.1). Then, for ej = Qhŭp — Up, there exists a constant C > 0 independent of h 
such that 


Illes ||| < chilislla- 
Proof. Taking v, = ej in error equation (4.3) yields 
aen, Cn) = a(QrÜn — Üh en) + li (ün, en) + Dg, en). 


According to Lemma 3.1, 3.4 and 4.4, we get 


les] < laCQndin — üns el + Is ins en)! + Hans en) 
< Bll] Qntin — às les] Chis] en] + Chill [es 
< Chllüs]l; es. || 


2 


which leads to the desired estimate. 


Let us derive the optimal order L? error estimate of eg by the standard duality argu- 
ment. The dual equation of (2.5) is seeking 9 € Hi AH? such that 


—P,(V - (AV)) = Pye; (4.10) 


where eg = Qoil, — ug as we defined in (4.2), and it follows from the elliptic regularity 
that [6], 


Illa < Clleoll. (4.11) 


Theorem 4.2. Assume that (H1) and (H2) hold. Let ti, € S? c H?N Hi be the solution 
of standard FEM in (2.4) and uj € V, be the solution of SFWG-FEM in (3.1). Then for 
eo = Qolly, — Up, there exists a constant C > 0 independent of h such that 


2 
lleoll  Ch^|[üg llo. 
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Proof. Testing Eq. (4.10) by eg yields 


lleoll? = M ^ (eoe): 


TES, 
= M (-V- V9), e0)7 
TER, 
= M (AVS, Veg); — 9, (AV®-n,e9 — ey) ar 
TEI, TEH, 
= 3 (AVS, Veg); — X (AVE -1, Qye0 — €5) ar 
TEX, TEX, 
= ps (AV® - n, eg — Quo) ar» 
TEA, 


where we utilize the information for the last second inequality that 


Y, (AV®-ny ep) 97 = 0. 
TES, 


Combining (3.1) with (4.5)-(4.7), we reach 


leol? = J, (AY, Vwenr— X, (A(V9 — Qu(V)) n, Queo — e) 


TEX, TEX, ee 
- M (AVG -n,e9 - Qseo)or (4.12) 
TEX, 


= a(ep, $) cm (9, en) E l)($, en) 
Let us bound the term |a(e;, ©)| first, 


a(e,, 9) = a(Qjüj — Up, P) + a(t, — Up, Qu) + aap, — un, 9 — QR?) 


4.13 
= a(QrÜn — Up, Ë) + li (ŭn, Qu) + D (ü5, Qn) + aC, — un. — Q9). ne 


In the following, we analyze these four terms in (4.13) successively. For the first term, 
we find the intermediate term (V, (Qj — tip), AQo(V®)), = 0, which comes from the 
Definition 3.2, the projection Qo and Green's formula 
(V, (Qui — tin), AQo( V 9));. 

= (V(Qoitn — tip), AQo(V9)) r + (2 ((Qptin — tin) — (Qottn — ün)), AQo(VS) - n) 3s 

= (V(Qolin — tin), AQo( V9)) r + (Quii — Qotin, AQo( V9) -n)55 

= - (Qoi — ün), V - (AQoC v) ),, + ((Qoün — tin) + (Qu — Qoi) AQoC V9) n.) .. 

= (Qui; — ty, AQo( V9) n) ar = O. 
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Applying the above equation, Lemma 3.2, 3.3, 3.4 and Cauchy-Schwarz inequality yields 


|a(Qpitn — tp, ®)| 


= » (AV, (Qiii — ün), Vw?) r 


TEX, 
< > | (AV, (Qui — tin), UC V9)) r| 

TE, 
= Y (Qi, — ün), AVE) r| 

f (4.14) 
RM | (v. (Quà, — tip), A(VO — Qo(V9))) / 

TEX, 

1/2 1/2 
«c p IV (Quis — wat (X [yes TOD 
TH, TE%, 


< Chl||Qnit, — ün |||llll2 
< Ch?llūnll2ll®ll2- 


For the second term, by means of Cauchy-Schwarz inequality, we have 


Ii Cà, Qn®)| 
=|} (A(vü, - Qvi) n. Qs(Q99) - Q38) 
TEX, 
< Y | (A(vain - Gà) n,Q- 9)),...| 
TES, 
SCS) [Vin - QV, lorllQo® — Fllar 
TEX, 
1/2 1/2 
ed Dy Il Vin, — 3 (x llQo$ — 1j 
TEH, TEH, 


Based on Lemma 3.2 combined with the trace inequality (4.8) and the inverse inequality 
(4.9), it is evident that 


Y [Vein — Qld, < chllüll2, 
TE%, 


Š Qo? - eid; < Chal. 
TEI, 


These two estimates give 
lli (45, Qn®)| < Ch? ll; |l» |l llo. (4.15) 
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For the third term, utilizing the projection operator Q,_; and Cauchy-Schwarz inequality, 
we get 


IL; (s, Qn®)| 


= »5 (AVü, - n,Qo$ — Q5(Qo9))5; 


TEX, 
< Y (Av, - O,-1( Vai) n, Qo® — Q4(Q99)),... 

T€4Z, 
€ C X Vin — Q, i (Vüg)llarllQo? — Q(Qo9)llor 

TEX, 
1/2 1/2 
«c (X llvà, — oaair] bs llQo$ — ast 
TEH, TEH, 


From the trace inequality (4.8) and inverse inequality (4.9) , it follows that 


2 Vän = QW AIMS 7 < Chllsll 
TES, 


Š Qo? - Q,(Q59)I2, < Chla. 


TEI, 
These two estimates lead to 
IL; Cds, Qj )| < Ch” là; llla. (4.16) 
For the fourth term, according to Lemma 3.1, 3.4 and Theorem 4.1, we have 


la(ü; — uy, 9 — Q1) 
< Bla, — us [ll — 2,8] 


z x (4.17) 
< CQàs = Qrin ll + Ileal- Ille -erll 
< Ch? |[i; |l; ||. 
Combined the estimates (4.14)-(4.17) together, we arrive at 
aes, &)| < Ch" [lyllallello. (4.18) 
Now, applying Lemma 4.4 and Theorem 4.1, we have 
Il (9, en)! < Chllella |]es|| < Ch? las llsllello, (4.19) 
Il, en) < ChIIell» |es]| < Ch? lapllallFlle- (4.20) 


Combing (4.12), (4.11) with (4.18)-(4.20), we obtain 


2 2 2un- 
lleoll^ < Ch^làs ||;]| ll; Ch |a; la lleoll; 
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which leads to the desired estimate. 


In the end of this section, we obtain the optimal order discrete H! and L? error estimate 
for the SFWG approximation (3.1). 


Theorem 4.3. Assume that (H1)-(H3) hold true and the exact solution u e H'^(Q) (0 < 
s € 1). Let up € Vj, be the solution of SFWG-FEM in (3.1). Then there exists a constant C > 0 
independent of h such that 

[fe — ual] © cre ll 


Proof. According to (4.1), Lemma 2.1, 3.4 and Theorem 4.1, we obtain 
[fu — ell < (lfe — às [] is — nal] | Qus — us] 
<C||Vu — Vügl| + Chl[ig l|; + Chllis l|» 
<Ch'||f ll. 


The proof is now completed. 


Theorem 4.4. Assume that (H1)-(H3) hold true and the exact solution u e H*^(Q) (0 < 
s € 1). Let uy = (ug, uy] € V, be the stabilizer free weak Galerkin finite element solution of 
(3.1). Then there exists a constant C > 0 independent of h such that 


llu — ugl| < Ch'^*|f |l. 
Proof. According to (4.1), Lemma 2.1, 3.2 and Theorem 4.2, we obtain 


llu — uoll S llu — à; l| + lün — Qottall + lIQods — uoll 
€ Ch |f || + Ch? lg la + CR? slo 
<Ch™ |f |l. 


This completes the proof. 


5. Numerical Experiments 


The goal of this section is to numerically verify the H!** convergence theory (0 < s < 1) 
for the stabilizer free weak Galerkin finite element method through three numerical exam- 
ples. In particular, we consider elliptic equations with various lower regularity and illus- 
trate the corresponding convergence rate. Let Q = [0,1] x [0, 1], and we use uniform trian- 
gular meshes J, where h denotes the spatial mesh size. Then the (P,(T), Po(e), [P5(T)]?) 
finite element space is 


V, := {v = fvg, v] : Vo € (T), vp € Pole), Vv € [P3(T)]?, VT € FH, Ve c OT]. 


The first example is presented to show the standard convergence rate for H?-regularity 
problem. The second and the third are supplied for studying the reliability of the SFWG 
method for H!** (0 < s < 1) problems with various regularity. Gaussian quadratures are 
used in the numerical integration. 
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Example 5.1. Consider (2.1) with the diffusion coefficient matrix 
x? +y? +1 0 
A= 2 2 š 
0 x*+y*+i1 
We choose the source term f and the boundary value g such that the analytic solution is 
u = sin(zx)sin(my ). 


Fig.1 shows the profile of the exact solution on the left and the stabilizer free weak 
Galerkin approximation on the right. We can find that the scheme is accurate. Table 1 
shows that the SFWG-FEM errors and the convergence rates of Example 5.1 in L?-norm 
and ||-||-norm for u. The convergence rates are of order O(h”) and @(h) respectively 
because the elliptic regularity is H?, which is consistent with the standard SFWG-FEM. 


Figure 1: The profile of exact solution (left) and SFWG-FEM solution (right) of Example 5.1. 


Table 1: Errors and convergence rates of Example 5.1 


h Wen — ull convergence rate ||Q,u—u,|| convergence rate 
1/8 3.6124E-01 — 7.0094E-03 — 
1/16  1.8130E-01 0.9946 1.7664E-03 1.9884 
1/32 . 9.0695E-02 0.9993 4.4252E-04 1.9970 
1/64 4.5344E-02 1.0001 1.1069E-04 1.9992 
1/128  2.2669E-02 1.0002 2.7676E-05 1.9998 


Example 5.2. Consider the problem (2.1) where the solution lies in the space H!**. Let 
the diffusion coefficient matrix A — I where I is a 2 x 2 identity matrix. We select the 
boundary value g and the source term f so that the analytical solution is 


u —x(1—-x)y( - yYK v x? -y2) 27, 0 «y x 1. 


For any0<e<y,lets=y—e>0. Then, ue Hin H1 5,0 «s <1. 
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We use the uniform triangular grids as before and carry out numerical computation 
with different step-sizes. For this example, a 32 x 32 mesh is used for plotting Fig.2. 
The profile of exact solution is presented in Fig.2 (left) and the solution calculated by the 
SFWG-FEM is shown in Fig.2 (right). We can find that the singular point is (0,0). 


Figure 2: The profile of exact solution (left) and SFWG-FEM solution (right) of Example 5.2 as y — 0.5. 


Since the solution u e H' nH ^ (9, as proved in the two theorems, we have the @(h*) 
convergence in ||-|||-norm and @(h'**) convergence in L?-norm, which are verified by our 
computational results. Table 2-4 show that the errors and convergence rates of Example 
5.2 as y = 0.25,y = 0.5 and y = 0.75. We can find that the H!-convergence rate is about 
of O(h9?5), o(h95) and O(Rh97), with L?-convergence rate being O(h!?°), O(h!) and 
O(h^), respectively. 

The numerical results are consistent with the theory for these cases and demonstrate 
that the SFWG-FEM method is also accurate and robust when the elliptic regularity is H!** 
(0 <s € 1), and a higher regularity leads to a faster convergence rate. 


Table 2: Errors and convergence rates of Example 5.2 with y — 0.25 


h Illu — uj convergence rate ||Q,u—u,|| convergence rate 
1/8 3.1102E-01 — 8.5605E-03 n 
1/16 2.6202E-01 0.2473 3.6025E-03 1.2487 
1/32 2.2025E-01 0.2505 1.5102E-03 1.2542 
1/64 . 1.8508E-01 0.2509 6.3335E-04 1,2537 
1/128  1.5556E-01 0.2507 2.6586E-04 1.2523 


Example 5.3. Consider (2.1) with the diffusion coefficient matrix 
"m x?- y?41 xy 
T xy x? +y? +1] 


We select the boundary value g and the source term f so that the analytical solution fits in 
with the one of the Example 5.2, which is 


u=x(1—x)y(1— y)(V x? + y2) ?**, 0<y<1. 


For any 0 < e < y, lets = y — e > 0. Then, ue Hin H!^,0 «s <1. 


16 Y. Wang, Y. Zou, X. Liu and C. Zhou 


Table 3: Errors and convergence rates of Example 5.2 with y — 0.5 


h ikem — ull convergence rate ||Q,u—u,|| convergence rate 
1/8 1.3165E-01 — 3.5534E-03 = 
1/16 9.4204E-02 0.4828 1.3075E-03 1.4424 
1/32 . 6.6960E-02 0.4925 4.7195E-04 1.4701 
1/64 4.7469E-02 0.4963 1.6873E-04 1.4839 
1/128  3.3612E-02 0.4979 6.0019E-05 1.4912 


Table 4: Errors and convergence rates of Example 5.2 with y — 0.75 


h Wenn — ull convergence rate ||Q,u—u,|| convergence rate 
1/8 6.8196E-02 — 2.0189E-03 — 
1/16 4.3163E-02 0.6599 6.8783E-04 1.5535 
1/32 .2.6779E-02 0.6887 2.2254E-04 1.6279 
1/64 1.6409E-02 0.7066 6.9931E-05 1.6701 
1/128  9.9693E-03 0.7189 2.1579E-05 1.6963 


We choose a symmetric positive definite diffusion coefficient matrix A and the solution 
lies in the space H!^(Y79, Fig.3 (left) and Fig.3 (right) show the convergence order of 
SFWG-FEM as u € H!°~¢ with 0 < e < 0.5 and u € H1 7?^* with 0 < e < 0.75, respectively. 
The blue dash line represents the convergence of ||-||-norm, the red dash line is of L?-norm. 
The red line in Fig.3 (left) is of slope 1.5 to be the reference of O(h!>~*) and the blue line 
is of slope 0.5 to be the reference of (12^ *). The lines in Fig.3 (right) are similar. 

We get the error and the convergence rate with respect to the mesh-size in Table 5-6, 
which also conforms well the theoretical analysis. This numerical experiment shows that 
the SFWG-FEM is still flexible and effective though the solution has low regularity. 


0 The convergence order of error " The convergence order of error 
10 i 10 ; 
E --4 
10" 2-g----üu7^ 4 
ee o OTe 
-2 
„10 
E 4 
10? 
- ar = 
E EL. di = [us — ull 
10 dere = |lu.—ull 5 
slope=0.5 
slope=1.5 
10? —— m 
10° 10 
Mesh size h Mesh size h 


Figure 3: The convergence order of SFWG-FEM as y = 0.5 (left) and y = 0.75 (right). 
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Table 5: Errors and convergence rates of Example 5.3 with y =0.5 


h ikem — up ||| convergence rate ||Q,u—u,|| convergence rate 
1/8 13210E-01 zT 3.3673E-03 — 
1/16 9.4264E-02 0.4869 1.2620E-03 1.4159 
1/32 6.6955E-02 0.4935 4.6173E-04 1.4506 
1/64 4.7462E-02 0.4964 1.6660E-04 1.4706 
1/128  3.3608E-02 0.4979 5.9594E-05 1.4831 


Table 6: Errors and convergence rates of Example 5.3 with y — 0.75 


h [nu — TAI convergence rate ||Q,u—u,|| convergence rate 
1/8 6.9643E-02 — 1.8109E-03 — 
1/16 4.3626E-02 0.6748 6.3529E-04 1.5112 
1/32 . 2.6943E-02 0.6952 2.1055E-04 1.5932 
1/64 1.6473E-02 0.7099 6.7315E-05 1.6452 
1/128  9.9948E-03 0.7208 2.1021E-05 1.6791 


6. Conclusion 


In this paper, we have proposed and proved the error estimates in both discrete Ht- 
norm and standard L?-norm for the stabilizer free weak Galerkin (SFWG) finite element 
methods for the second order elliptic problems with low regularity solutions. This analysis 
helps to remove the compulsory H?-smoothness assumption for the real solution in all the 
related literature. We assume that the solution of elliptic problems has H ! 5-regularity, and 
we have proved that the H!-convergence rate and the L?-convergence rate are of s-order 
and (1+s)-order, respectively. The (P(T), P._1(e), [Pes 7 i) elements are employed for 
SFWG-FEM. Our strategy is taking H? finite elements as the intermediate approximation. 
We approximate the elliptic problem using the traditional FEM with at least H? smooth 
bases, then approximate this standard finite element solution via SFWG-FEM. Therefore, 
the (P,(T), Po(e), [P2(T)]*) elements are reliable to construct weak Galerkin finite element 
bases for the low-regularity problems. We also give three numerical examples for this. 

Our result is an important supplementary for the new stabilizer free weak Galerkin 
finite element method theory. We believe this result can be extended to more types of 
PDEs that the finite element method can be applied to and more general meshes. 
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